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Abstract 

We study a free interface problem related to combustion of condensed matter and 
some non-equilibrium exothermal phase transitions. In spite of a variety of non-trivial 
dynamical scenarios exhibited by the model the solutions are uniformly bounded and 
the interface velocity is a smooth function. The main result of the paper establishes 
existence of a compact connected attractor for the classical solutions of the problem. 
Numerical evidence leads to the conjecture that the fractal dimension of the attractor 
is finite. 

Submitted to Interfaces and Free Boundaries 

1 Introduction 

This paper presents a study of attractors for a two-phase Stefan problem with kinetics. We 
show that classical solutions of the problem approach a compact connected attractor in the 
uniform norm. We also demonstrate via direct numerical simulations that the attractor has 
a finite correlation dimension. 

The free-boundary problem that is the subject of the paper arises naturally as a math- 
ematical model of a variety of exothermic phase transition type processes, such as solid 
combustion \W\ also known as Self- propagating High-temperature Synthesis or SHS fK)\ . 
solidification with undercooling ^^1, laser induced evaporation ^3], rapid crystallization in 
thin films |231 etc. These processes are characterized by production of heat at the interface, 
and their dynamics is determined by the feedback mechanism between the heat release due 
to the kinetics and the heat dissipation by the medium. In addition to its theoretical in- 
terest SHS has industrial applications as a method of synthesizing certain technologically 
advanced materials for high-temperature semiconductors, nuclear safety devices, fuel cells 
etc., see [20], and also PU| for a popular exposition. SHS propagates through mixtures 
of fine elemental reactant powders (e.g., Ti + C, Ti +2B), resulting in the synthesis of 
compounds. 
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There is a substantial literature that treats analytical aspects of the initial-boundary 
value problem for different sharp-interface models with kinetics related to the problem 
()2. 1112 .411 below, see [iHl 123 [211 13 EHl IHj- These works are concerned with basic issues 
of mostly local in time existence. See also ^3 IIH] where a finite interval version of the 
problem with linear kinetics and the asymptotic behavior when a kinetic parameter tends 
to zero are discussed. Dynamics on an semi-infinite spatial interrval for a one-phase model 
are investigated in |^. The principal result of this paper is asymptotics for the position 
of the front s{t) of the form s{t) ~ kt" + 0{t^) for t ^ oo, < /3 < a, where a = 1/2 
or 1, depending on the value of undercooling. Needless to say that these asymptotics, 
being important on their own right, cannot capture order one variations in the temperature 
profile and velocity. These variations, their compact structure and presumably low Hausdorff 
dimension are the subject of our work (the graph of the time history of solutions in Fig. 1 
below gives some idea about their complexity). We also note recent papers by Brauner et 
al., pP-j2]! which study dynamical behavior of solutions of a related problem. In particular 
they consider perturbations of traveling-wave initial data and investigate their instability 
and bifurcations. 

The objective of our work is to investigate asymptotic behavior of a propagating front. 
A variety of complex asymptotic dynamics (cf. ^) arise only in the context of a infinite 
spatial interval. This necessarily calls for a problem on an infinite spatial domain, note 
that for a bounded domain the flame front reaches the external boundary in finite time 
and extinguishes. It is important to realize that the image of a ball under the evolution 
is not compact for any finite time. The potential theory computations in Sec. 3-4 allow 
us to extract the part of the evolution that compactifies (the contribution from the free 
interface), while the heat losses force the contribution from the initial conditions decay with 
time exponentially. We believe that this clear structure of the attractor indicates that careful 
potential theory estimates are unavoidable if one is to study asymptotic patterns and prove 
existence of a compact attractor in our case. 

This study is partially motivated by numerical simulations and in particular by the nu- 
merical experiments described in [7] , where it was demonstrated that the system generates a 
remarkable variety of complex thermokinetic oscillations. The dynamical patterns exhibited 
by the system, as the governing parameters are varied, include a Hopf bifurcation (see the 
rigorous proof ^^), period doubling cascades leading to chaotic pulsations, a Shilnikov-Hopf 
bifurcation etc. These patterns are well-known for the finite-dimensional dynamical systems 
and hint at the possibility that the essential dynamics of the free-interface problem may be 
finite-dimensional as well. 

At the same time, in jH] we demonstrated that a 3 x 3 system of ODEs obtained as a 
pseudo-spectral approximation to the one-phase free-boundary problem exhibits dynamics 
that mimics that of the infinite-dimensional system to a surprising degree. For the one-phase 
problem we were able to prove that compactness and finite dimensionality of the attractor 
do take place [TT] . 

It should be noted that the two-phase problem is somewhat more physically sound than 
its one-phase counterpart and appears in various applications. However, the methods of 
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the papers dealing with the one-phase problem are not directly applicable to the two-phase 
Stefan problem with kinetics which is the subject of the present communication. This is due 
to the additional temperature field behind the propagating interface (in the product phase) 
being not easily controllable. We overcome this difficulty in the present paper and show that 
the temperature behind the interface is sufficiently well-behaved to render compactness of 
the attractor. 

Having proved compactness one is naturally curious as to how "large" is the attractor 
in terms of some appropriate measure. Currently we are not able to answer this question 
analytically due to, as we believe, purely technical difficulties although we have little doubt 
that the dimension is finite for the two-phase case. As the measure of complexity of the 
attractor and therefore of the asymptotic regime, we use the correlation dimension introduced 
by Grassberger and Procaccia (see, for example, JH])- The correlation dimension is based on 
the idea that if the evolution can be described by a finite number of degrees of freedom then 
the time series of observations on the system should be spatially correlated. This spatial 
correlation is measured by the corellation sum, which is directly related to the integral of 
the standard pair correlation function c(r): 



If the dimension is v then it is easy to see that C(l) ~ for small I. We compute the 
correlation dimension of the attractors generated in the direct numerical simulation of the 
problem. 

The rest of the paper is organized as follows. In Sec. 2 we state the free-boundary prob- 
lem and present some minimal background information on local existence and uniqueness. 
Sec. 3 and 4 constitute an analytical core of the paper. Here we employ the single-layer 
representation for the solution to obtain a natural decomposition of the solution into two 
contributions, one from the initial conditions, and another one from the free boundary. We 
obtain certain potential theory based estimates for both contributions to the solution and 
its spatial derivative, which are instrumental for the proof of compactness in Sec. 5. The es- 
timates are proved to be uniform with respect to the sup norm of initial data, assuming that 
the nonlinearity in the kinetic free-boundary condition satisfies some natural requirements. 

In Sec. 5 we utilize the estimates that show that the contributions from the free boundary 
are uniformly bounded and decay at infinity. Together with the uniform bound on the 
spatial derivative, they allow us to apply a version of the Arzela-Ascoli theorem, which 
guarantees that the contributions from the free boundary for initial data from a fixed ball 
form a precompact set. We complete the proof of existence of a compact attractor by 
using an appropriate abstract result from dynamical systems. Finally in Sec. 6, from direct 
numerical simulations of the free-interface problem we estimate the correlation dimension of 
the attractor. 
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2 The free-boundary problem. Local existence and 
uniqueness. 

We study the following free boundary problem: find s(t) and u{x,t) such that 



Ut = u, 



'XX 



7M, X ^ s{t) 



(2.1) 



u{x, 0) = Uq{x) > 



(2.2) 



g[u{s{t),t)]=v{t) fort>0, 



(2.3) 




and the derivatives and u~ are taken from right side and left side of the free interface 
respectively. The last term in the heat equation ()2.H) is due to the heat losses into the 
medium surrounding the combustible or solidifying substance via Newton's cooling law with 
a non-dimensional coefficient 7 > 0. 

The surrounding matter is assumed to be at the temperature of the fresh combustible 
mixture at —00 (the original phase in the phase transition interpretation). By the same 
token the heat loss will reduce the temperature in the product phase to that of the medium. 
Thus the behavior of the solution at infinity should satisfy \\mx^±ao t) = 0. It should be 
remarked that the presence of the heat losses 7 > only improves the analytical properties 
of the solutions. For 7 = the boundary condition at 00 should be replaced by convergence 
to a constant. 

The dynamics of the physical system are determined by the feedback mechanism between 
the heat release due to the kinetics g[u\x=s{t)) and the heat dissipation by the medium. To 
illustrate the meaning of the kinetic term, for example, in the context of solidification, we 
note that for some substances in the presence of strong supercooling of the original phase the 
phase transition temperature measured at the interface may deviate considerably from the 
equilibrium one and is functionally related to the interface velocity. This dependence called 
the interface attachment kinetics can be different for different substances due to various 
microscopic mechanisms responsible for the incorporation of the product at the interface 
into the crystalline lattice. 

The second interface condition ()2.4|) (the Stefan boundary condition) expresses the bal- 
ance between the heat produced at the free boundary and its diffusion by the adjacent 
medium. As the problem describes propagation of the phase transition front, the first in- 
terface condition (j2.3|) is a manifestation of the nonequilihrium nature of the transition; 
its analog for the classical Stefan problem is just u\x=s{t) = 0. We should mention that in 
contrast with the nonequilibrium problem, the dynamics of the classical Stefan problem is 
relatively trivial. 
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To discuss properties of the attracting set for classical solutions of the free interface 
problem ()2.1|) - ()2.4|) we need to first establish their existence and uniqueness for all times. 
We say that u{x,t),v(t) form a classical solution of ()2.H1 - ()2.4|) if 

(i) u{x,t) and v{t) are continuous for t > 0; 

(ii) u^x and Ut are continuous for x 7^ s(t), t > 0; 

(iii) Equations ()2.1|) - ()2.4p are satisfied. 
We quote the result as it was stated in 6 : 

Theorem 1 Suppose that the kinetic functions g satisfies the following assumptions: 
{Al) g{u) is a continuously differentiable, monotone decreasing, negative function on (0, 00) 
with g{0) = — fo for some velocity —vq < 0; 
{A2) g{u) is suhlinear: Imiu^oogiu) /u = 0; 
and that the initial data uo{x) > are bounded. 

Then there exists one and only one classical solution u{x, t) > and v{t) of the free interface 
problem ll2.1\ )- \2.4\ )- The solution is uniformly bounded for all t > 0. 

For the reader's convenience we outline the scheme of the proof. First, the problem is 
reduced to an integral equation for the interface velocity using the single layer potential 
representation: 

G{x,t,C,0)uo{OdC- / G(x,t,s(r),r)e-^(*-")t;(r)rfr, (2.5) 
■00 Jo 

where 

G(x,t,e,r) =exp{-i^-^}[47r(t-rrV2 (2.6) 
4(t - r) 

is the heat kernel and s{t) = j^v{T)dT. 

Taking the limit of (|2.5p as x — ^ s{t) and using the kinetics condition (j2.3p . we obtain an 
integral equation in terms of v only: 



-7t 



G(s(t),t,e,0)no(Of^e 



t, s(r), r)e-^(*-")t;(r)dr, (2.7) 



Next, we show that for sufficiently small time intervals, the mapping K defined by the 
right hand side of the integral equation 

ft 



v{t) = g 



.-It 



G{s{t),t,^,0)uo{Od^ 



G{s{t),t,s{T),T)e-''^'-^^v{T)di 



Kv 



(2.8) 

defines a contraction on an appropriately chosen closed set of continuous functions which 
yields local existence. We also remark that the velocity v{t) can be shown to be an infinitely 
differentiable function. The centerpiece of the global existence proof is the a priori estimate, 
which allows us to extend the local solution indefinitely. 
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3 A priori estimates: spatial decay of solutions 

In order to demonstrate existence of an absorbing set we need to establish spatial decay 
of the interface contribution to the classical solutions of the problem. First we obtain an 
estimate for the solution on the interface. 

Theorem 2 Letu{x,t),v(t) be a classical solution of / fF7|) - (|i:^.^| j and \ \uo\ \ = sup_^^^^^\uo{x)\ < 
M, then 

\u{s{t),t)\ <2M + Rg, \v{t)\ < g{2M + Rg) (3.9) 
where Rg and A are constants dependent solely on the kinetic function g. 

Proof. First we prove the estimate for the interface temperature: iplt) = g{v{t)) = 
Ii — I2, where /i and I2 are the two parts of the right hand side of ()2.7j) . 
It is rather obvious that < Me~'^*: 

f°° e""^* (s(t) - f)^ 



|e-^* / G(s(t),t,e,0)no(Orfel = ;r^ll%ll / exp{-^^^— = e-^*||«o|| 

J-00 Zy/tlT J-00 4^ 

Now, since the kinetic function satisfies the condition (A2), for any e > there exists 
f 1 > such that \g{ip)/4'\ < £^ if gii') < —vi. We subdivide I2 as follows: 

G{s{t),t,s{T),T)e-''-'-^^g{^{T))dT= [ +[ = h + h, (3.10) 



where X\ = {^^1 - ^^i < fi'(^(^)) < -^^o, < r < t} and X2 = (0, t)\xi- 
For /a we have: 

G{s{t),t,s{T),T)g{4j{T))dT\ <Vi\ [ G{s{t),t,s{T),T)dT\ 

xi Jo 



00 



Here we have used the observation that \s{t) — s(r)| = |f (0I(^ — ''")> VQ{t — r) for some 

T <i<t. 

Now, let us interpret Ii{t) = Pipit) as a mapping; then the following estimate holds for 
its norm: 

IIP^II < I [ Gisit),t,siT),T)gii;iT))dT\ 

•JX2 

< I / G{s{t),t,s{T),T)etp{T)dT\ 

JX2 

< e\m\ r G{s{t),t,s{T),T)dT\<-U\\. (3.11) 

Jo 
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Since ^ = 1^-1^- P^, we have + PVll = \\h{t) - h{t)\\ < M + vi/vq. On the other 
hand, by choosing e = fo/2, we have ||P|| < |, and therefore 

\il){t)\ = u{s{t),t) < 2M + 2vi/vo. 

The constant 

Rg := 2vi/vo (3.12) 
is the constant referred to in the statement of the theorem. Simuhaneously, 

\v{t)\ < g{2M + Rg) (3.13) 

■ 

Now using the uniform bounds for the interface velocity and temperature that we have 
just estabhshed it is easy to obtain the uniform estimate for the entire field through the 
maximum principle: 

\u{x,t)\ <2M + Rg, (3.14) 

From now on we shall assume that g{u) is a monotonically decreasing differentiable 
function on [0, oo] with \g'\ < C (recall that the velocity v = g{u) is negative) and 

■^0 < ~g{'^) ^ K) for some Vq, vq > 0. (3.15) 

Both conditions are satisfied for the standard Arrhenius kinetics that in appropriate rescaled 
and normalized variables has the form ()6.3p . The existence of the lower bound Vq in particular 
seems to be crucial for the uniform boundedness of solutions. 

For the compacness result we need certain decay estimates for the contribution from the 
free boundary (see ()3.17|) - ()3.18|1 below) which are established in the the following 

Lemma 3 Let u{x,t),v{t) be a classical solution of l\2.1\ )- \2.4\ ) with \ \uo\ \ < M then for the 
contribution from the free interface 

^{x,t)= [ G{x,s{T),t-T)e-^^'-^^v{T)dT 







the following estimates hold: 



\^>{x,t)\<^ (3.16) 
2V7 



(a) For X < s{t) 



m..«)|<^A^exp{-Hit^_tz^} (3.17) 

V f +47 2 4t 
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(Hi) For X > s{t) 



Vq 

exp(— a(x — s(t)), for X — s(t) > 2Vo/ 'J 

V7 



mx,t)\<{ V/ (3.18) 

forO<x-sit)<2Vo/j 



where a = mm(t>o/4, 7/2Vb). 

Proof. First we obtain a very simple bound, which is vahd for any x, t: 

t (x._,(^))2/4(t_^) 

V47r(t - r) 

< Vo / e-^(*-) ^ , dT<Vo e-^'-' 







Note that this estimate is very different from the one obtained through the maximun prin- 
ciple, ()H.14j) : the dependence on the norm of the initial conditions is absent. 
Ahead of the interface x < s{t) 



t _ ^ ^_(^._s(r))V4(t~r) 

f (r)(ir| 



-7(t-r)e 



v^47r(t - r) 



< Vo / e-^^*-"^exp{-^ y ' y \ '''''' '-^} 



4(t - r) 

~ i/tt 2 4t 

exp{-^^-4 - 7 ^ - ^ } ; 

t;o|x-s(t)| \x-s{t)\'^ 

- / 2 , A ^^P^ 9 77 ^ 

'vl + 47 2 4t 



To estimate the free-interface contribution to the solution |\I^(a;,t)| behind the interface 
X > s(t) we split the interval of integration into two subsets: = {r G [0,t] : s(r) < 
+ x)/2} and its complement X2 = ^ [0^^] ■ ■s(t) > + x)/2}. 

J' G{x,t,s{T),T)e-^'^'-^^ \v{r)\dT = J + J = h + l2, 

Xl X2 
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For the first integral we have 

J 2^7r{t - r) ' ^ ^' 

Xi 

Xi 

- ^ " r)-V^exp[-(a; - .(t))|]e-(*--)cir 

Xi 

(a;-s(t))/(2i;o) 



erf exp[-(.-.W)^] < J% cxp[-(. - 



The following inequahties 



(t - r) ~ ' 2 ' (t - r) ~ ' 2 ' x - s(t) 

have been used to replace the exponent in the Gaussian kernel, which gave rise to the 
exponential decay factor. 

For the integral I2 we obtain 

r ^ -l{t-r) ^^_ {x-s{T)f 

J 2J7r(t-T) Mt-r) J' ^ ^' 



X2 

00 



J 2y/Tfr] yJ-K J 

{x-s{t))/{2Vo) ^{x-s{t))/{2Vo) 



exp(-7(a; - s{t)) / {2Vo)), for 7(0; - s{t))/{2Vo) > 1 

<{ 7o"' (3.19) 
for < 7(x - s{t))/{2Vo) < 1 

The final incqualties in the above estimate are based on the following elementary observa- 
tions: if aVb > 1 then 

/•OO 2 7*00 /»00 

/ exp{—brj^)dri = —= / exp{—rj^)dr] < —= / r] exp(—bri^}dr] = — ^exp(— 6a^); 
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on the other hand 



exp{—bri'^)dri < I exp(~bri'^)dri = ^ 

Jo 2Vb 



Thus for X > s{t) we obtain 



|^(x,t)| < < 



Vo 
2v^ 



^exp[-(x-s(t))j] + 



exp(-7(x - s(t))/(2K,)), for 7(x - s(t))/(2Vo) > 1 



Vo 
2^7 



for < 7(x - s(t))/(2K)) < 1 



(3.20) 



Obviously the direct contribution from the initial conditions is bounded by the norm of 
the initial conditions: 



-7* 



G{x,t,^,0)uo{Od^< 



-It 



2Vt7r 



exp{-^^^K = e-^*|KII (3.21) 



4 Estimate for the derivative 

The proof of compactness is based on a version of Arcela-Ascoli theorem and uses an estimate 
for the derivative of the solution. Via differentiation of the representation of the solution 
fl2.5|l . the derivative for x ^ s{t) is expressed as follows: 

u^ix, t) = -e-^' / G^ix, t, 0)uoiOd^ + / G^ix, t, s(r), r)e-^(*-^)^;(r)cir, (4.22) 



Lemma 4 Let v(t) be a continuous function on [0,T], define the derivative of the boundary 
contribution as 



$(x,t) = ^ [ G(x,s(T),t-T)e-^^'-^^v(T)dT 
ox Jo 

Then for every < t < T |$(x, t)| < const 



(4.23) 



Proof. The estimate ahead of the front, i.e. for x < s{t), is treated as follows. In the 
estimates ahead of the interface we replace exp(— 7(t — r)) by 1. Consider separately two 
cases: \s(t) — x\ > 1 and \s(t) — x| < 1. 
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For the case \s{t) — x\ > 1 



X 



s(t) e-(^-<-^))V4(t-r) 



< 4/e 

^^^^ ^ (x - .(t))^ + 2(x - - s{t)) + (.(t) - .(r))% v{T)dT 



8(^-t) V^(*-^) 

< C U _ ^-)-3/2g-(a^-s{t))V8(t-r)g-«ok-.(t)|/4g-^2(i_^)/8^^ 

^y^^-v,\.-s{tm j-t _ |. ._3/2 -(.-.(t))V8(t-r)^ 

32\/2\/oe-^ol^-*WI/^ /"^ _„2 , lev^^o e-"ol^-'(*)l/^ 

< 1= — ^ — / e ' dr] < -— — (4.24) 

ey/7rvQ\s{t) - x\ Jo cvq \s{t) - x\ 

In the last estimate we used the following simple observations: ^e"^ < 1/e, for ^ = 

- — -. — > 0, \s(t) - x\ > \s(t) - x\. \s(t) - x\ > \s(t) - s(t)\ > vo\t - t\ and sub- 
8{t-r) 

stitution T] = \s{t) — x\(t — t)^^^'^/^/8 to obtain the error function integral. 
For the less involved case — a;| < 1 we split the integral into two parts 

ft X- s(t) e-(^-*W)'/4(t-T) 

< r k - s{t)\ + \s{t) - s{t)\ e-(^-^W)V4(t-r) 
~ ' Jo 2(i - r) ^in{t - t) ^' 

< / M!l_iK^_lL_e-(-K*))V4(t-r)^^ + / ^_--dr 



A Jo 4 ^V^Jo y/{t - r) 

<Yl + JLr e-^di = ^(1 + ^) (4.25) 



2 Vo^/tt Jo 2 tiQ 



where = (a; - s{t)){t - t)-^/^/2 and ^ = Voy/{t-T)/2. 
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Thus, collecting estimates in ()4.24|) - ()4.25|) we observe that the derivative ahead of the 
interface x < s{t) decays exponentially 



|$(a;,t)|< (4.26) 

^(1 + ^), s{t)-l<x<s{t) 
2 vq 

The part of the estimate concerning the derivative behind the interface x > s{t) causes 
some difficulties due to the fact that x can be close or even equal to s(r); it is treated as 
follows. For our purposes it suffices to prove the estimate for t starting from a certain t > 0. 
Assume for convenience that t > 1. For the second integral in ()4.22|1 we split the interval of 
integration into two subsets: Xi = {'^^ ^ [O?^] • — x\ > \s{t) — x\/2} and its compliment 
X2 = [0,t] \ Xi- Note again that s(t) is a monotone function. We have 

Gi:{x,t,s{T),T)e-^^'-^^v{T)di 



Xi X2 

For Ii we use two subsets of Xv Xi = Xii U X12 where x > s(r) for Xn x < s(r) for Xu- 
For the integrals we get respectively 

{x-s{T)f 

/■|x-s(r)| """"P^ 4(t-r) ,,,,, 
Til = / — X ^ ^ — e ' X \v{Tj\dT 



<l \G^{x,t,s{T),r)e-^'^'-^^\x\v{T)\dT (4.27) 
= /i + /2, (4.28) 



Xii 

00 



2(t-r) 2y^ 

00 

- ^/ l^-^WI(^-^)"'^'^^p[-(^-^W)'i^(rr7)]^^ 


00 

< [e-^'d^ = Vo (4.29) 
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(note that on Xu \x — s{t)\/2 < \x — s{t)\ < \x — s{t)\). In the above integral we have made 
a substitution \x — s{t)\{t — r)~^/^/4 = ^. For the second part 



'12 



\X — S{T) 

2(t - r) 



exp[ 



(x - sir))' 



X 



4(t-r 



Xl2 



{x - s{r)r 



)_e-l{t-r) X \y^r)\dT 



< 



< 



\s{t)-s{r)\ 
2(t-r) 



exp[ 



4(t-r 



Xl2 



^e-^(*^") X |t;(r)|c/r 



r 



— ^(ir 



1/2 



(4.30) 



We remark that a time independent estimate for the above integral is also valid. 
For I2 we shall replace a function of the type xe~^ by its maximum l/v^2e 



\x — s r 



exp[ 



2{t 



X 



- 4r) f 
A{t - r) 



X2 

< 



2^Tx{t - r) 



X \v{r)\d'] 



x-s{t)\ {x-s{T)f 

X exp[— ^-T- — ^4^1 X 



X2 

1 



2v/(t^ 
dr 



4(t-r) ' 2y^{t-r) 



dr 



'ire 



< 



-^0 



dr 



< 



^0 



(t-r)- ' J \sit) - x\/2 - V0V2 



ne 



X2 



X2 



(4.31) 



We have also replaced (t — r) in the denominator by its minimum \s{t) — x\/2 \\v\\ on X2^^d 
observe that meas{x2) — l^i^) —x\/vo. 

By collecting the estimates in ()4.29p - ()4.3ip we obtain a uniform bound for the derivative 
behind the interface x > s{t), 



\^{x,t)\<Vo + i^+ ° 



4^7 vq^/2t^ 

This concludes the proof of the lemma. ■ 

We also note that the first term in ()4.22|) can be easily estimated, 



(4.32) 



G5(x,t,e,0)Mo(Oe-^*rfe 



< e 



-7t 



exp[ 



X 



-7* foo 



< Ikol 



2^1 



X exp[- 



At 



2Vvrt 



l^o(OMe 



-7* 



At 



2Vt 



< \\uo\ 
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since the integral on the last line is equal to 2 riexp{—r]'^)dr] = 1, with r] = \x — ^\/2^/i. 
Thus we get the following 

Corollary 5 For any t > to, to > 0, the derivative of the solution is uniformly bounded: 



Ux{x,t)\ = 




(4.33) 



We reiterate that the proof of compactness in the next section is based on a version of 
Arcela-Ascoli theorem and uses the estimate for the derivative of the free-interface contri- 
bution, < const. Thus the estimates of this section are crucial for the compactness 
result. 



5 Absorbing set and attractor 

In this section we use the estimates obtained above to establish existence of a bounded 
absorbing set and of the attractor which is compact in the space of continuous functions. 
It can be easily verified that all the estimates and analytical properties of the solutions can 
be obtained without the heat losses. On the other hand the problem with the heat losses 
exhibits uniform exponential decay in time of the contribution from the initial data which is 
utilized in the proof of compactness of the attractor. 

The integral representation ()2.5p describes the evolution of the initial temperature dis- 
tribution uq: u{t) = T{t)uQ. We think of the evolution as taking place for the functions on 
(— oo, oo) in the moving coordinate system attached to the free boundary x' = x — s{t). Note 
that all the results in Sec. 3-4 are in tems of x'. 

It is convenient to split the semigroup operator T into two parts: the contribution of the 
free boundary 

Ti{t)uo{x') = - [ e~"'^'~^^G{x' + s{t),t,s{T),T)v{T)dT (5.1) 
Jo 

and that of the initial data 

/oo 
Gix' + sit),t,^,0)uoiOd^ (5.2) 
■oo 

We rephrase the esimates in ()3.21|) and ()3.16|) as the following result that establishes 
existence of an absorbing set for the evolution. 

Proposition 6 (i) The semigroup T2 is uniformly exponentially contracting in C : 

sup ||T2(t)M°|| <exp(-7t)A^ 

for any ball 

X = {ueC; ||u||<A^} 
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(ii) For any e > 0, the ball Babs := {u E C : \u\ < Vo/{2^) + e} is an absorbing set for 
all bounded subsets of C for the semigroup T. Here the radius of the absorbing ball reflects 
the contribution of the free interface alone. 

Next we prove that the boundary contribution to the evolution, i.e. the operators Ti{t) 
are uniformly compact. Namely, the following proposition holds: 

Proposition 7 For any to > the orbit of the ball 0{X,to) = Ut>toTi{t)X is relatively 
compact in C. 

Proof. For the version of Arzela-Ascoli theorem appropriate for C it is sufficient to 
have uniform boundedness for the derivative and uniform decay of the family of functions 
as \x'\ — > oo. From Lemm^we see that for the initial data in a ball, the contribution from 
the interface exhibits a spatial decay uniformly with respect to time. On the other hand, 
the estimate (??) demonstrate that the spatial derivative is uniformly bounded. 

Now it is a simple matter to construct a finite e-net for 0(X, to)- First we choose a finite 
interval ~L < x' < L, beyond which the functions of the family are smaller than e, it is 
possible to accomplish because of the uniform in time spatial decay. In view of the uniform 
bound on the derivative, the functional family is equicontinuous. Therefore the restricion of 
0(X, to) on to [—L,L] is compact by the regular Arzela-Ascoli theorem. By extending the 
elements of the e-net from [—L, L] to the whole line by zero we obtain an e-net in 0{X, to)- 
■ 

The properties of the evolution operator T(t) described in the above propositions allow 
us to apply the abstract general result (see, for example, |23l Chap. 1]) that in our situation 
can be stated as follows: 

Theorem 8 The uj-limit set A of the absorbing set Ba is a global exponential compact at- 
tractor for the metric space C ; A is the maximal attractor in C and it is connected. 

6 Numerical dimension of attractor 

Having proved compactness one is naturally curious as to how "large" is the attractor in 
terms of some appropriate measure? For the one-phase problem we have been able to prove 
that the Hausdorff dimension of the attractor is finite ^T]. However, due to the additional 
temperature field behind the propagating interface (in the product phase), methods of the 
papers dealing with the one-phase problem are not directly applicable to the two-phase prob- 
lem. Currently we are not able to overcome these, we believe, purely technical difficulties, 
although we have little doubt that the dimension is finite for the two-phase well. 

The question arises also, whether the presence of the temperature field behind the front 
affects the "size" of the attractor in comparison with the one-phase case (see ^2])- We 
provide an answer to this question by computing the correlational dimension of the attractors 
generated via direct numerical simulation of the problem. 
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While the Hausdorff dimension is convenient for analytical estimates, it is highly non- 
trivial to compute and requires too much storage and CPU time. More convenient com- 
putationally is the correlation dimension. Although in general dcorr < dnausdorff, they are 
usually very close. We follow the now standard procedure for computation of the correla- 
tion dimension [16j. Namely, consider the set {f/j, i = 1,...A^} of points on the attractor 
Ui = U{T + ir), where T ^ 1. We consider a discreet approximation of the solution in 
the space M'^ by sampling the solution at k points Ui = {u{xi,T + ir), ...,u{xk,T + ir)). 
We measure the spatial correlation between the points on the discreet approximation of the 
attractor with the correlation integral 

1 

C{1) = lim —-{number of pairs with the distance p{Ui, Uj) < 1} 

If for small Z, C{1) scales as l'^ then the correlation exponent v can be taken as the correlation 
dimension of the attractor dcorr- For practical calculations the frequency of sampling r, the 
number k of points in space where the solution is sampled at each time, and the number of 
samples N are determined empirically. Similarly, for the low sample dimension k a better 
approximation for dcorr may be obtained if the Euclidean distance p is modified by inclusion 
of a weight. 

To obtain a numerical approximation of the attractors we solve the initial value problem 
()2.1|) - ()2.4|) for sufficiently large time until the asymptotic regime is attained. Obviously 
the dimension of the attractor should not depend on the choice of initial data, which was 
confirmed by direct numerical simulations. Problem 1)2.11) - ()2.4|1 was solved in the frame 
attached to the free boundary on a finite interval [— i^, L\ with the Dirichlet conditions 
w{~L, t) = simulating the decay of the solution at — oo, and dw{L, t)/dx = corresponding 
to the stabilization of the temperature in the product phase. According to our observations 
the results are practically insensitive to the increase in the interval length after L ~ 30 (see 
[Zj for the details of the numerical algorithm). We remark that in contrast with the one- 
phase case where L ~ 10 was sufficient, one needs a rather large spatial domain to obtain 
consistent numerical results. 

To represent different dynamical regimes we use the Arrhenius kinetics , 

V = g{u) := -exp[a — ^ ], (6.3) 
a -\- [I — a)u 

where (in the context of combustion) a is proportional to the activation energy (Zeldovich 
number), and a is the temperature ratio of the fresh mixture and the product, see e.g. |Z]. 

Thus, the attractor is represented as a set in M*', where k is the number of sampling 
points of temperature profiles. We choose time snapshots of the solution for every 0.08 in 
the interval of the asymptotic regime (200 < t < 1800) and consider them as a discrete 
approximation of the attractor in W^. The correlation dimension for this discrete set is 
evaluated as explained above. As a control experiment we selected a periodic asymptotic 
regime, a = 4.5. It is immediately confirmed that dcorr ~ 1 as one should expect. 

In contrast, for a = 5, a = 0.05 the regime is chaotic as is illustrated in Fig. 1 that 
presents a series of snapshots of spatial temperature profiles. One can see from the log-log 
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Figure 1: Time history <t < 4000 for chaotic dynamics u{x, t) vs. x, t. 



graph of the correlation integral (Fig. 2) that in this case dcorr ~ 2. From our observations 
on a variety of regimes it appears that the dimension cannot be much higher than 2. 

7 Concluding remarks 

The compactness result has been proved here in the presence of heat losses for any nonzero 
heat loss. Although we chose to operate in spaces of continuous uniformly bounded functions 
on the infinite interval, we are convinced that compactness can be established even for zero 
heat loss if spaces with weaker topology are used, for instance in the space of continuous 
functions bounded on each finite interval. 

Results of this paper have been proved for the kinetic function satisfying the bounds 
in (??). These bounds are quite physical and cover a wide range of important applica- 
tions. Nonetheless, our numerical experimentation with different types of kinetic functions, 
including unbounded ones demonstrate that the asymptotic dynamics are insensitive to the 
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Figure 2: Correlation integral for 7-, 8- and 9-point samples. 



behavior of the kinetic function for large temperatures. On the other hand, our results from 
[Hj provide global existence and uniform boundedness of solutions for a wider class of kinetic 
functions, namely for sublinear kinetics. Therefore we strongly believe that the principal 
result of this paper holds for this CclSG clS well. 

It is interesting to compare the proof of the compactness above to that for the one-phase 
problem JUj. Although the estimates in the two-phase case are more involved due to the 
presence of the temperature field behind the propagating interface, once they are obtained, 
the representation of the evolution semigroup is more transparent than in the one-phase 
case. 

Also, it is rather remarkable that the numerical estimates for the correlation dimension 
of the attractor above and for the one-phase case (see [12]) yield roughly the same value. 
Indeed, such an outcome is rather unexpected because the two-phase problem seems to 
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possess more "degrees of freedom" than its one-phase counterpart. 

Finally, for the one-phase problem we have been able to prove that the Hausdorff dimen- 
sion of the attractor is finite [llj. Currently we are not able to overcome certain difficulties 
that we believe are of purely technical nature but it appears safe to conjecture that the 
dimension of the attractor is finite for the two-phase problem as well. 
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